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Abstract 

We consider the massive tricritical Ising model Al(4, 5) perturbed by the thermal oper- 
ator (^1^3 in a cylindrical geometry and apply integrable boundary conditions, labelled 
by the Kac labels (r, s), that are natural off-critical perturbations of known confor- 
mal boundary conditions. We derive massive thermodynamic Bethe ansatz (TBA) 
equations for all excitations by solving, in the continuum scaling limit, the TBA func- 
tional equation satisfied by the double-row transfer matrices of the lattice model 
of Andrews, Baxter and Forrester (ABF) in Regime III. The complete classification 
of excitations, in terms of (m, n) systems, is precisely the same as at the conformal 
tricritical point. Our methods also apply on a torus but we first consider (r, s) bound- 
aries on the cylinder because the classification of states is simply related to fermionic 
representations of single Virasoro characters Xr,s[Q)- We study the TBA equations an- 
alytically and numerically to determine the conformal UV and free particle IR spectra 
and the connecting massive flows. The TBA equations in Regime IV and massless RG 
flows are studied in Part II. 

1 Introduction 



Ever since their introduction ^ |^ , Thermodynamic Bethe Ansatz (TBA) equations have 
been an important tool in the study of both massive and massless integrable quantum field 
theories. Although extensive studies have been carried out on scaling energies of vacuum 
or ground states only relatively few excited states 0, ||, 0, § have proven amenable to 
TBA analysis and these are primarily restricted to massive and diagonal scattering theories. 
So despite considerable successes, the application of TBA methods has been hampered by 
inherent limitations. The primary obstacle is that to date there is no systematic and unified 
derivation of excited state TBA equations. Indeed, the current treatments of excited states 
are at best ad hoc and fall well short of a complete analysis of all excitations. Here we address 
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these limitations and propose a systematic approach based on the lattice. More specifically, 
we will show in a series of papers that both massive and massless excited TBA equations 
can be systematically obtained by studying the continuum scaling limit of the associated 
integrable lattice models. Perhaps the most important input from the lattice approach is 
an insight into the analytic structure of the excited state solutions of the TBA equations. 
Previously this structure had to be guessed. In stark contrast, in the lattice approach, the 
analyticity structure can be probed by direct numerical calculations on finite size transfer 
matrices. 

Although the methods developed in this paper are very general, for simplicity and con- 
creteness, we consider as a first example the massive tricritical Ising model (4, 5) perturbed 
by the thermal operator (fi^s- Although this is a non-diagonal scattering theory and more 
complicated from the viewpoint of integrable quantum field theory, this is the simplest case 
beyond the Ising model and Lee- Yang theory for analysis by the lattice approach. The in- 
tegrable lattice model associated to the thermally perturbed tricritical Ising model is the 
interacting hard square model or generalized hard hexagon model solved by Baxter H, IT^ 



This model, with its Z2 sublattice symmetry, is known to be in the universality class of the 
tricritical Ising model. More generally, this model is the special case L = 4 of the Al lattice 
models of Andrews, Baxter and Forrester with the L = 3 model being the usual Ising 
model. Generically these Al models, with their Z2 height reversal symmetry, are in the 
universality class of multicritical Ising models. Moreover, the continuum scaling limit of the 



Al models realize the <y9i 3 thermal perturbation of the si{2) unitary minimal models [|T3 
and the ground states are described by the Al-2 TBA equations of Zamolodchikov 0. 

There have been many relevant studies of the tricritical hard square or A4 lattice model 
and the more general Al models from the lattice viewpoint. For the model, the off-critical 
TBA functional equation for periodic boundary conditions has been derived and solved [|1^, 
r5[ for the bulk properties and correlation lengths. The off-critical TBA functional equations 



for the Al models were derived by Kliimper and Pearce |]T6|, 0, |T8[ . But only the critical or 
"conformal TBA" equations were derived and solved in the critical scaling limit for the central 
charges and conformal weights. The very same off-critical TBA functional equations for Al 



models were subsequently derived |19| in the presence of integrable boundaries showing 
that the TBA functional equations are universal in the sense that they are independent 
of the boundary conditions. A biproduct of introducing boundaries is that the problem of 
classifying the excitations becomes much easier. This is reflected in the fact that at criticality 
the cylinder partition functions are given as linear forms in characters rather than the usual 
sesquilinear form on the torus. Indeed, by a judicious choice of (r, s) boundary conditions, 
the cylinder partition function is just a single Virasoro character Xr,s{fl) and the complete 



classiflcation of excitations [20| in terms of (m, n) systems [21, 22| is related to a fermionic 



representation of the character Xr,sic[)- This simpliflcation enabled pO| the complete analytic 
calculation of the conformal cylinder partition functions of the A^ model with 6 different 
conformal boundary conditions (r, s) which are conjugate to the 6 primary flelds of the 
tricritical Ising conformal fleld theory. The generalization of these results to the critical Al 
lattice models with (r, s) boundaries is currently in progress . 

In 1998, Pearce and Nienhuis |2^ analysed the off-critical continuum scaling limit of the 
TBA functional equations for the Al lattice models to derive from the lattice the full set 
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of massive and massless Al_2 ground state TBA equations conjectured by Zamolodchikov. 
The scale parameter mR is simply related to the scaling limit of lattice parameters by 



fi = ^= hm Nf (1.1) 



or more precisely 



Af^oo,t^O 



At" 

R = lim Na, m = lim — (1-2) 



where a is the lattice spacing, m is a mass, R is the continuum length scale, = (L + l)/4 is 
the correlation length exponent and t = is the deviation from critical temperature variable 
with p the elliptic nome appearing in the Boltzmann weights of the Ai models. 

Our primary goal in this series of papers is to extend the analysis of Pearce and Nienhuis 
to all excitations of the A4 model both in the massive and massless regimes. This entails 



perturbing the analysis of O'Brien, Pearce and Warnaar pO| off criticality. To handle the 
problem of classification of all the excitations it is easier to introduce boundaries and to 
work on the cylinder even though a study of boundary properties is not our primary goal. 
Our immediate goal is to study the flow of excitation energies from the UV {R = 0) to the 
IR [R — > 00) limit. In the massive case we are thus able to compile the conformal-massive 
dictionary that eluded Melzer In the massless case considered in paper II [^], we follow 
the renormalization group flow from the tricritical to the critical Ising model fixed points. 
This leads to a flow between the characters of these theories. 

In this paper we consider just the massive regime. The layout of the paper is as follows. 
In Section 2 we define the A4 lattice model. We then describe the classification of excita- 
tions and discuss their unique labelling in terms of quantum numbers. This classification 



is in fact identical to the classification ||20| at the tricritical point. In Section 3 we present 
the derivation of the off-critical massive TBA equations. The numerical solution of these 
equations is presented in Section 4. Throughout we concentrate on the three (r, s) boundary 
conditions with s = 1 rather than presenting exhaustive results for the six distinct (r, s) 
boundary conditions. We believe these results are indicative of what can be achieved. While 
the calculations are similar for the other boundary conditions, we point out that in some 
cases there are further subtleties related to the appearance of frozen zeros. We conclude 
with a general discussion. 



2 Lattice Approach 
2.1 lattice model 

The A4 RSOS lattice model is defined on a square lattice with spins or heights a = 1, 2, 3, 4 
restricted so that nearest neighbour heights differ by ±1. This model corresponds to the 
special case L = 4 of the A^ RSOS models of Andrews, Baxter and Forrester |Tl| with spins 
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L and Boltzmann weights (A = 7r/5 for the model) 
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Here u is the spectral parameter, A = 7r/(L+l) is the crossing parameter and - 
is one of the standard elliptic theta functions as given in Gradshteyn and Ryzhik 



2p^l* sin M ~ cos 2u + p'") (1 - p'") 

n=l 

00 
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n=l 



Integrability derives from the fact that these local face weights satisfy the Yang-Baxter 
equation. 

The elliptic nome p is a temperature-like variable. The Ai lattice models are critical 
for p = and off-critical for p 7^ so we introduce the deviation from critical temperature 
variable t = p^. There are four distinct off-critical physical regimes depending on the sign of 
u and t: 



Regime I: 
Regime II: 
Regime III: 
Regime IV: 



-7r/2 + A < n < 0, 
-7r/2 + A < M < 0, 
< M < A, 
< M < A, 



-1 < t < 
< t < 1 
< t < 1 

-1 < t < 0. 



(2.^ 



It is convenient to express the nome p in terms of a real parameter e > by 



p 



le 



Regimes II and III 
Regimes I and IV 



so that IpI < 1 and 

t = p^ = ± exp(— 27r£:). 

In particular, the elliptic 'di functions satisfy the quasiperiodicity properties 

§i{u + 7r,p) = -■&i{u,p) 
•di{u — ilogp,p) = — p^"'^e~^*"7?i(M,p). 



(2.9) 

(2.10) 

(2.11) 
(2.12) 
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Regimes III and IV are of interest in this series of papers since they are associated, in the 
continuum scaling limit, with the massive and massless thermal perturbations of the unitary 
minimal models respectively. In Regime III, considered in this paper, p is real whereas in 
Regime IV p is pure imaginary. Regimes I and II relate to I^l-i parafermions and so are not 
considered here. 

Although we will not use it in this paper, an alternative formulation of the A4 model 
is the particle or T2 tadpole representation as shown in Figure 1. This formulation is ob- 
tained by folding the A4 diagram and identifying the states related by the Z2 height reversal 
symmetry. More specifically, we can identify the states a = 1,4 with /i = 1 and regard this 
as indicating the presence of a particle or an occupied site and we can identify the states 
a = 2, 3 with /i = and regard this as indicating the absence of a particle or a vacant site. 
Once we fix the sublattice of the square lattice which has odd heights, the identification of 
height and particle states is a one-to-one correspondence. The adjacency constraint on the 
heights of the model translates into the exclusion of simultaneous occupancy of adjacent 
sites by particles. In this way the particle representation is seen to be equivalent to a model 
of interacting hard squares on the square lattice. The Boltzmann weights of this hard square 
model are simply given by replacing the heights a = 1,2,3,4 by the corresponding particle 
occupation numbers = 0, 1. 



Figure 1: A^ and T2 graphs of allowed neighboring states for the lattice model and its 
particle (hard square) representation. 

The isotropic phase diagram of the A^ or interacting hard square model is shown in 
Figure 2 alongside the corresponding renormalization group flow in the continuum scaling 
limit about the tricritical point. The continuum scaling limits in Regimes III and IV give 
rise to the massive and massless flows respectively where the perturbation parameter is mR 
or /i. 

2.2 Double row transfer matrices 

To ensure integrability of the A4 model in the presence of a boundary [|1^ we need commuting 
double row transfer matrices and triangle boundary weights which satisfy the boundary Yang- 
Baxter equation. The triangle weights for the (r, s) boundary condition with 1 < r < 3 and 
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Figure 2: Isotropic phase diagram of the or interacting hard square model showing the 
phase boundary (sohd and thick sohd curves) between the fluid and square ordered solid 
phases. The fugacity of the particles is z and L is the attraction between particles. The 
tricritical point T separates the line of Ising critical points (solid curve) from the first-order 
line (thick solid curve, Regime III). The model is exactly solvable on the first-order line 
and its analytic continuation (dashed line, Regime IV). The exact solution curve and Ising 
critical line are tangential at T. The renormalization group flow in the continuum scaling 
limit about the tricritical point T is shown below for comparison. The central charge c is 
indicated at the tricritical, critical and trivial fixed points. 



1 < s < 4 are given by 



K 



( ^ r ± 1 


= ( 







?9i(rA) ) ^\{\) 



(2.13) 



m) 



The parameters ^l, are arbitrary and can be taken to be complex, however here we restrict 
them to the real interval (A/2, A). To obtain conformal boundary conditions at the isotropic 
tricritical point u = A/2, t = p = we should choose |2^ = = A/2. Integrability in 
the presence of these boundaries derives from the fact that these boundary triangle weights 
satisfy the boundary Yang-Baxter equation. 

The face and triangle boundary weights are used to construct [ITO] a family of commuting 
double row transfer matrices D{u). For a lattice of width iV, the entries of D{u) are given 
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diagrammatically by 
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(2.14) 
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It is further convenient to define the normahzed transfer matrix 



where 



and 



with 



t{u) = Sr^s{u)S{u) 



S{u) 



Diu) 



SrJu) 



hr{u) 

hs{u) 



^?i(2m-3A,p)?9i(2m + A,p) 

-iyhr{u - iL)h-r{u + ^Ljhsiu - iR)h-s{u + ^/j) 

^9i(A,p)^9i(m + (3 - r)A,p)i9i(M + (1 - r)A,p) 

— A,p)'i9i(M + 2A,p) 
di{\p)d^{u + (3 - g)A,p)^4(M + (1 - s)\p) 
'di{u,p)'di{u — X,p)'d/i{u + 2A,p) 



(2.15) 

(2.16) 
(2.17) 

(2.18) 
(2.19) 



It can then be shown that the normahzed transfer matrix satisfies the universal TBA 
functional equation 



t{u)t{u + X) = / + *(« + 3A) 



(2.20) 



independent of the boundary condition (r, s). This is precisely the same TBA functional 
equation that holds in the periodic case |l^. Since the transfer matrices commute this 
functional equations also holds for each eigenvalue t{u) of t{u). 

The TBA functional equations will be solved for the finite-size corrections to the eigen- 
values D{u) of the double row transfer matrices D{u). In the scaling limit, the finite-size 
corrections to the eigenvalues D(u) are related to the excitation energies E{R) of the asso- 
ciated perturbed conformal field theory by 



-\ log D{u) = Nf^^M + Ksiu) + E{R) + 0(1; 



f2.211 



where /buik(w) is the bulk free energy, br,s{u) is the boundary free energy and the anisotropy 
angle is given by 



(L + l)u, Regimes III and IV 

1^=( 2(L + l)u ^ , 

— , Regimes I and II. 

±j — 1 



(2.22) 



7 



Depending on the boundary conditions, the boundary free energy br^s{u) may contain an 
interfacial free energy contribution. The bulk and boundary free energies can be calcu- 



lated p9|, ^ ^ by the inversion relation method. Despite the appearance of correc- 
tions, the system is not in general conformally invariant. The system is conformal however 
at critical points which can occur in the ultraviolet {R 0) and infrared {R — >■ oo) limits 
with 

^^^-^ + Ar,. + n, neN (2.23) 

where c is the central charge of the appropriate conformal field theory. A,, ^ are the related 
conformal weights and n = 0, 1, . . . labels the tower of descendants. The largest eigenvalue 
occurs for the vacuum or ground state with the boundary condition (r, s) = (1, 1). In this 
case Ai 1 = and n = 0. The massive R oo scaling limit in Regime III, however, is trivial 
in the sense that for this ground state RE{R) — » corresponding to c = and the scattering 
of free massive particles. 

2.3 Classification: {m,n) systems and quantum numbers 



TBA functional equations admit infinite families of solutions for the eigenvalues t{u). The 
analyticity properties are therefore crucial in selecting out the required solutions. The trans- 
fer matrix eigenvalues D{u) are entire functions of u and are characterised (up to an overall 
constant) by their zeros in the complex u plane. It is precisely at these zeros that logt(u) 
is non-analytic but analyticity of logt{u) is required to solve the TBA functional equations 
by Fourier series. From quasiperiodicity, the matrix t{u) and eigenvalues t{u) are doubly 
periodic. It follows that the eigenvalues t{u) are doubly periodic meromorphic functions. It 
is convenient to fix the period rectangles as 

panod rectangle =/ ("i |) ><(-?. R^^"- (2.24) 

I (— |, ^) X {—7Ti6,TTiej, Regime IV 

so we then only need to consider the analyticity inside these period rectangles. In Regime IV 
there is an additional symmetry within the period rectangle 

t{u±'K/2 + 'Kie) =t{u) (2.25) 

so we can restrict ourselves further to the rectangle (— f, 2A) x {—nie, nie). The normalization 
factors relating D{u) to t{u) only introduce extra zeros and poles on the real axis. Since 
D{u) is real symmetric for real m, that is D{u) = D{u)'^, it follows that, for any eigenvalue 
D{u) or t{u), the distribution of zeros in the upper and lower half planes are identical and 
simply related by complex conjugation. It is therefore sufficient to classify the eigenvalues 
by the patterns of zeros in the upper half period rectangle. 

It turns out that in Regime III the pattern of zeros inside the periodic rectangle is 
qualitatively the same as in the critical case. This was observed by direct numerical diago- 
nalization of a sequence of finite-size transfer matrices approaching the scaling limit N oo, 
t —>■ for modest sizes of A^. As a consequence, in Regime III we can use the known clas- 
sification pO] of eigenvalues at the tricritical point [R = 0) in terms of (m, n) systems and 
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this classification will apply for any R in the range < i? < oo. This simplifying feature 
does not hold in the massless Regime IV where in stark contrast some of the patterns of 



zeros qualitatively change under the flow as R increases [^. Since the classification of ex- 
citations in Regime III of interest here is the same as at the tricritical point we summarise 
the salient features of this classification here. We limit discussion to the boundary condition 
(r, s) = (1,1). The other cases are similar |^ although in some cases it is necessary to 
introduce two (m, n) systems for a given (r, s) boundary condition. 

The analyticity properties of D{u) are relevant in the two analyticity strips 

A 3A 

< Re{u) <—, 2\< Re{u) < 4A (2.26) 

We refer to these as strip 1 and 2 respectively and label them by i = 1,2. Since the 
Boltzmann weights are real and positive for < m < A, strip 1 is referred to as the physical 
analyticity strip. From direct numerical diagonalization of D{u) with the (1, 1) boundary 
condition we observe that, apart from a pair of zeros on the real axis at m = A + ^ and 
M = 5A — ^ induced by the left boundary triangle weight, each eigenvalue has zeros on the 
lines Re(M) = —A/2, A/2, 3A/2, 2A, 4A corresponding to the edges and center lines of the 
two analyticity strips. Specifically, 1-strings and 2-strings occur in strip 1 and 2. A 1-string 
is given by a single zero Uj in the center of a strip such that 

A 2-string is a pair of zeros {uj, u'j) on the edge of a strip with equal imaginary part and 

(Re(.,),Re(.;))4^-'/''''/')' (2.28) 
^ ^ ^ \(2A,4A), strip 2. ^ ' 

Distributions of zeros for two typical eigenvalues of D{u) for = 16 are depicted in Figures ^ 
and We note that for finite A^ the 2-strings do not fall precisely on the lines given in 
equation (|2.28|) , but that this deviation decreases exponentially, that is as \S\^ with \S\ < 1, 



as A^ increases. These patterns are consistent with the crossing D(u) = D(X — u) and 
transpose D{u) = D{u)'^ symmetry of the double row transfer matrix. 

Given an eigenvalue, we denote the number of strings in the upper half period rectangle 
as follows: 

rrii = number of 1-strings in strip i = 1,2 , ^ 

^ ^ ' (2.29) 

Tii = number of 2-strings in strip i = 1,2. 



The relations |^0|] between these numbers determining the string content take the form of 



an (m, n)-system ||2T|, |22 



m + n = ^{Nei+Im) (2.30) 

where m = {mi, 7712), n = {ni,n2), ei = (1,0), and X is the A2 incidence matrix with 
entries Xj^k = S\j-k\,i- Clearly here we require that mi,m2 and A^ are even. For the leading 
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AX 



Figure 3: Zeros within a period rectangle in the complex li- plane of the largest eigenvalue 
of D{u) with string content m\ = 6, ni = 3, m-i — 2,n2 — 1. 



21 



3X 



4X 



Figure 4: Zeros within a period rectangle in the complex u-plane of another typical eigenvalue 
of D(u) with string content mi = 6. ni = 3, ni2 = 2,n2 — 1. This pattern of zeros is obtained 
by permuting the ordering of the 1-strings and 2-strings within each strip of Figure 3. 
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excitations mi,m2,n2 are finite but ni ~ N/2 as N ^ oo. Indeed, the vacuum or ground 
state is given by mi = m2 = n2 = and ni = N/2. 

For each system size N, there are many eigenvalues with the same string content (m, n). 
These eigenvalues are distinguished by the relative vertical orderings of the 1 and 2-strings 
within the period rectangle along each strip. Denoting the imaginary parts of the 1-strings 
in strip i by < f f ^ < ■ ■ ■ < Vm- , and the imaginary parts of the 2-strings in strip i by 
< < ■ ■ < Wn} , we see that in Figure |^ 



< ) < v^^^ < V? 
whereas in Figure ^ which has the same string content, 

0<t;f)<^f)<t;f. 



(2.31) 



(2.32) 



a) (i) 

Notice that the 1-strings Vj and 2-strings Wj labelled by j = 1 are closest to the real 
axis. Clearly, the total number of possible orderings for a given string content (m, n) is 
i^rn^^) {"^^m^^^) " Summing over all allowed string contents (m, n) using (|2.30|) then gives 

= (2.33) 



(m.,n) 



y-v f mi + ni\ f 1712 + n2 



nil J \ m.2 



This is indeed the correct number of eigenvalues as given by the dimension of the double 
row transfer matrix, which is A^^, the number of A^-step paths from 1 to 1, where A is the 
A^^ adjacency matrix. 
In the scaling limit 

^ = !!^= lim Nt"^ (2.34) 

4 Af-^oo, t^O 

the positions of the 1- and 2-strings grow logarithmically as 



v^//\w^-^ ~ —-logt + const = -logA^ log yU + const. (2.35) 

J J 5 5 



More specifically, we define the scaled locations of the strings in strips i = 1, 2 as 
I3f = hm Uvf + 7 logtV J = 1, 2, . . . , mi 



" 7vJ^t-.o ( + ^ ) , = 1, 2, . . . , ma 



(2.36) 



An excitation with string content (m, n) is uniquely labelled by a set of quantum 
numbers 

/ = (/(^) |/(^)) = (If \ /«,..., 4\) I If ), If , . . . , /S) (2.37) 
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where the integers J-*^ G N with i 



w 



1, 2 give the number of 2-strings whose imaginary parts 



^ are greater than that of the given 1-string Vj . Clearly, the quantum numbers satisfy 



n,>/P>4^^>...>/^)>0 



rO') 



1,2. 



(2.38) 



Conversely, given the quantum numbers we can read off the values of mi and 1712 and then 
ni and n2 are uniquely determined by the {m,n) system. For given string content (m, n), 
the lowest excitation occurs when all of the 1-strings are further out from the real axis than 
all of the 2-strings. In this case all of the quantum numbers vanish /j*'' = 0. Bringing the 
location of a 1-string closer to the real axis by interchanging the location of the 1-string 
with a 2-string increments its quantum number by one unit and increases the energy. At 
the tricritical point where R = these jumps in energy are quantised in a tower with energy 
levels given by |20] 



E 



lim 



RE{R) 



27r 



c 1 ^ 

1 — mCm 

24 4 



2 mi 
\^ \^ 



i=\ j=l 



(2.39) 



Here C is the A2 Cartan matrix and the central charge is c = 7/10. 

More succinctly, the generating function for the finite-size spectra of a cylindrical lattice 
built by M applications of the double row transfer matrix with (1,1) boundary condition is 
given by a finitized Virasoro character [pTl B2 



TTTiCm 



mi + rii 




m2 + n2 


mi 




m2 



E [m,n) 

where the sum is over the finite (m, n) system, q is the modular parameter 

q = e-^^w^i'^^ 



(2.40) 



(2.41) 



and M/N is the aspect ratio of the lattice. In the isotropic case when u = X/2, the anistropy 
angle d = ti/2 and the geometric factor sind = 1. The g-binomial or Gaussian polynomial 
is defined by 

m+n 



m + n 
m 



1+...+/,, 



/i =0/2=0 



{(l)m{q)n 

0, 



m, n > 



(2.42) 



otherwise 



with the g-factorials {q)m = (I""?) ■ ■ ■ (1~?™) for rn > 1 and (g)o = 1. In the limit q —>■ 1 the 
g-binomials reduce to the usual binomial coefficients and the partition function just counts 
the number of states as in (^1331) . Note also that 

(2.43) 



lim 


m + n 




00 






n— >oo 


m 




m 



After using the {m,n) system to eliminate ni and n2, the finitized character gives the 
fermionic representation of the usual Virasoro character in the limit N 00 



lim xZ\(l) 



-c/24 



qjmCm 



mi , r?i2 even 



mi 



\mi 
m2 



Xi,i(?)- 



(2.44) 
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More generally, as explained in EO], the finitized partition functions for (r, s) boundary 
conditions at i? = are given in terms of finitized characters 



(m,n) 



n 

1=1,2 



m + Hi 

rrii 



XSH?) (2.45) 



where the (m, n) system for the (r, s) boundary condition is 



m + n = -{Nei + Xm + Ss-i + e4_ 



(2.46) 



Hence in the limit ^ oo we recover the known conformal cylinder partition functions in 
terms of Virasoro characters 



(2.47) 



3 TBA Equations in Regime III 

3.1 Sector (r, s) = (1, 1) 

In this section we derive the TBA equations for the (r, s) = (1, 1) boundary by solving the 
TBA functional equations in the scaling limit for even A^. We follow closely the derivations 
in 



20|| and ||2J]. The derivation for other boundary conditions is similar. We begin by 



factorizing t{u) for large A^ as 



tiu) = f{u)g{u)l{u) 



(3.1) 



where f{u) accounts for the bulk order- A^ behaviour, g{u) the order-1 boundary contributions 
and l{u) is the order- 1/A^ finite-size correction. We will solve for /(m), g{u) and then l{u) 
sequentially. 

For the order- A^ behaviour the second term on the RHS of the TBA functional equation 
(|2.2(]| ) can be neglected giving the inversion relation 



f{u)f{u + X) = l. 



(3.2) 



The prefactor in ( 2.15|) induces poles of order 2A^ at m = 2A and m = 4A, and a zero of order 
2A^ at M = 3A. The required solution with this analyticity is 



f{n) 



L'7?i(f ,t5/4) 



-| < Re(M) < f 



2N 



3A 



9A 



2 < Re(u) < 2 



(3.3) 



Similarly to the critical case, putting this solution into the TBA functional equations 
implies the order-1 functional equations for g{u) 



g{u)g{u + \) 



-I < Re(n) < f 



l + g{u-2\), 



3A 



9A 



2 < Re(M) < 2 • 



(3.4) 
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To solve for g{u) we need to take into account the order-1 zeros and poles introduced by the 
order- 1 pref actor in (|2.15|) 

The order-1 zeros of D{u) cancel exactly the poles of Si^i{u). However S{u) introduces poles 
at M = -| + i^p, f + 2A + i^, 4A + i^, and double zeros at m = f + i^p, 3A + 
where p = 0, ±1. Thus the solution in strip 1 is given by 



9{u) 



^i(5(m- A/2)/2,t5/8) 
i?2(5(M- A/2)/2,t5/8) 



(3.6) 



The solution for g{u) in strip 2 is more involved and requires solving the functional 
relation for log g{u) by Fourier series. To proceed we fix lines of constant real part in the 
centers of each of the two strips in the w-plane and a real coordinate x by 

.= + ^' ^^"P^ (3.7) 
[3A + f, strip 2. ^ ^ 

It is then natural to define generically for the functions h = t, f,g,l the notations 

h,{x) = h{^ + ^), |Im(x)|<7r (3.8a) 

tx 

h^ix) = h(3X + |Im(x)|<7r (3.8b) 
5 

Hi{x) =1 + hi{x), H2{x) = 1 + /i2(x). (3.8c) 
In the variable x, the functional relations become 

9i{x - —)gi{x+ —) = 1 (3.9a) 

71% TCI 

g2ix--)g2{x + -) = Giix). (3.9b) 

One can show that the ratio g2{x) / gi{x) is free of zeros and poles for |Im(a;)| < vr. Similiarly, 
Gi{x) is analytic and non-zero in |Im(x)| < ^. Thus solving the functional relation for g2{x) 
using Fourier series, we find 

log5f2(a;) = logfi'i(2;) + e* logG'i(x) (3.10) 
where the kernel in the convolution is 

= (3.11) 

and z/ = 5/4. We do not need the explicit solution g2{x) since we only need to evaluate it in 
the scaling limit. 
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The functional relations for the finite-size corrections l{u) are obtained from ( |2.2(J| ) using 
(|33) and O 

h{x - -)h(x + y) = Tt(x) (3.12a) 
«--?)M-H-y) = ^. (3.12b) 

To solve for log/i(a;) and log/2(a^) we need to remove the singularities arising from the zeros 
in the interior of strips 1 and 2, that is, the mi and 1112 1-strings {^^iv^p} and {3A±if^^^}. 
Using elementary solutions of the inversion relation l[x — ^)l{x + y) = 1 with a single zero 
inside the analyticity strip we find 

y ^9l(f + |^t;f ,t5/4) ^^(1 _ |^^),t5/4)' 

is free of zeros and poles inside each strip. Applying Fourier series to the logarithms of 
( |3.12| ) and using Fourier inversion thus gives the following nonlinear integral equations valid 



for |Im(x)| < TT 



log/,(x) = ^log ^ J. ]l I +e*\ogUx) + C, (3.14a) 

log h{x) = 5J log I . 5, ^ ' I (2) 5. + ^ * log ^ + (3.14b) 

where the integration constants Ci and C2 are multiples of tt? related to the choices of 
branches for the logarithms. By going to the critical limit t — > 0, and comparing these 
equations with the corresponding equations in [^], one finds that Ci = C2 = 0. In effect we 
are fixing the branches of the logarithms of Ij exactly as in the critical case. 
We next write these equations in terms of ti{x) and t2{x) as 

log ti (x) = log /i (x) + log gi {x)+e* log T2 (x) 



™i (ix I 5„- ,(1) „Q fix 5„- ,(1) 



5 > 



+ 5^1oj 



(3.15) 



log t2{x) = log /a (x) + log 5(2 (x) + e * log Ti (x) - e * log Gi{x) + 
^y^^^ ^i(f + |^t^f,t!)^,(f -|^t;f,tt) 
fc°^^2(f + |^4'^^^)^2(f->f,t!)- 

We are interested in the solutions of these equations in the scaling limit. Replacing tt 
by jj: we see that all dependence on t disappears and only a dependence on remains. We 
assume the relevant functions have the general scaling form 

/i(x) = lim /i(x + logA^) (3.16) 

A''— +00 
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and set 



= lim U{d - log^) = - log/i), 



1,2 



(3.17) 



with /i = 2Hi > 0. The ei('(9) are precisely the pseudo-energies and "i? is the rapidity. 

Now taking the scahng hmit of the nonhnear integral equations using (p.3|), ( p.9| ) ( |3.14| ) 
and (|3.17|) gives the excited TBA equations 



logtanh - - 2^1og^tanh (2 + 2 ^ ^^2 ~ 2 ^ 

i=i 

, log(l + e-^2{^')) 



mi 



(1) 



1 

2^ 



cosh(^9 — 



9 "^2 Q n{2) Q 0(2) _^ 

e2(i9) = 2mi?cosh^9 - logtanh^ - - ^log[tanh (- + tanh (- - 

fc=i 



(3.18) 



0/ log(l + 



-ei{f')^ 



cosh('i9 — -d') 



where we have moved to the scaled locations of zeros /Sj*"* as given by (2] 



The excited TBA equations contain mi + m2 extra parameters which are the locations 
of the zeros inside strips 1 and 2. These can be determined by considering the scaling limit 
of the TBA functional equations 



TT TT 
ti{x - l-)ti{x + 1-) = 1+ t2[x) 

h{x - t-)t2{x + = 1 + ti{x). 



(3.19) 
(3.20) 



Setting X = Y + Sfj*'' we see that the LHS must vanish. Hence one can show that in the 



scaling limit 



t2(/?f^-y-log/^)- , 



(1) ■ 
-1 = e"^ 



j = 1,2,... ,mi 



(2) 



(2), 



-log/i) = -1 = e'^fc = 1,2 



(3.21) 
(3.22) 



where n^p and n^^^ are odd integers. Moreover these integers, which are determined by 
windings, must be precisely the same as in the critical case {R = 0), namely. 



n 



n 



(1) 
j 

(2) 



2(mi - j) - ma + 1 + 2/, 



(1) 



2(m2 - A;) - mi + 1 + 2/^- 



(2) 



j = 1,2, . . . ,mi 

= 1,2, . . . ,1712. 



(3.23) 
(3.24) 



Applying ( |3.17| ), the auxiliary conditions determining the locations of zeros become 

e2{(3f -^)=nf\t, j = l,2,... ,mi (3.25) 
=4'^^^, A; = 1,2,... ,m2. (3.26) 
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For numerical purposes we need a more explicit form of the auxiliary equations obtained 
by replacing with (3j — y in the TBA equations 



2m/2sinhM^^ = / — 



log(l + e-^^(^9)) 



7712 



2ir siiih(/3'" - tf) 



En 

k=l 



tanh f ^ + ^JL^l^ 
V 4 2 



k=\ 



V 4 2 



(2) , ,^(1) 



■ 1 I 1 2 / ^' 

I log tanh — 



(1) 



V 4 



rft? log(l + e-^2(;^9)) 



z ^ log 



mi 

i^log 



2vr sinh(/?i'^ - ^9) ^.^^ 
- i log 



4 2 



(3.27) 



tanh(^-^^' 



(2) , /^(l) 



tanh^ - ^ 



(2) 



V 4 



We propose ( |3.18| ), together with the auxiliary equations ( p.27| ), as the TBA equations for 
all excitations in the massive perturbation of the tricritical Ising model on a cylinder with 
the (1, 1) boundary condition. 

It remains to relate the finite-size corrections to the pseudo-energies to obtain the scaled 
energies. Again following ||2^ , one can determine the finite-size corrections to the eigenvalues 
of the double row transfer matrix from ( |3.14| ) as 



1 
2 



mR cosh X 



N 



-^logDi(x) = ^2cosh/3j^^ - ^ I rft9cosh?91og(l + e-^2('')^ 



2^ 



(3.28) 



where we have neglected terms of order o(l/A^). The scaling energies of excitations are 
therefore 



RE{R) = 2mR cosh /?] 



(1) 



mR 
~27 



/oo 
di!} cosh log (1 + e" 
-oo 



(3.29) 



3.2 Analysis of UV and IR limits 

One can check that the UV limit i? — of the massive TBA equations reproduces the 



"conformal TBA" equations of O'Brien, Pearce and Warnaar |20]. To do this, one should 
make the identifications 



•d ~ log 



and for the pseudo-energies 



mR ^(i) , mR a) 



, , mR 

e{x) ~ e I log — — h X 



(3.30) 



(3.31) 



Doing this we find that this limit indeed reproduces the known "conformal TBA" equations. 
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The IR limit i? — » oo of the massive TBA equations can also give many insights on 
the field theoretic behaviours. In this limit one can interpret the auxiliary equations (|3.27|) 
as Bethe ansatz equations for mi massive particles and m2 massless particles interacting 
with each other. The momenta and energies of the particles determined by the equations 
will introduce corrections in the total energy corresponding to vacuum polarization due to 
a large but finite value of mR. 

In this limit we find that the TBA equations become 



ei(^) 



2mi?cosh?9 + 0(l) 



log 



cosh?? 



mi 



cosh -i? + 1 1 , 



(1). 



cosh + cosh p, 
cosh — cosh 13 ^ 



and the first auxiliary equation is simplified as 



(3.32) 
(3.33) 



2mR sinh /3 



(1) 



R' 



This implies 



(1) 



[m2 + 1 + 71-^)71 

2mR 



Substituting these results into the second auxiliary equation in ( p. 27]) we find 



{nil + 1) log 



-i sinh /3f ^ + 1 
—i sinh — 1 



(2) . 



Solving this in the IR limit gives the limiting locations of zeros in strip 2 



sinh 13 



(2) 



cot 



(2) 



A; = 1,2, 



.2(mi + 1) 

Finally, the large R limit of the scaling energy is given by 

m r°° 

E{R)r^2mim / dd cosli^ S{'d)e 

where S{'d) is given by 0(1) term in 62 ("(9) 



,m2. 



2mR cosh i? 



S{^) = tanh^ - 



m2 

n 

k=l 



(2). 
k 



cosh 'd — cosh (3^ 

(2) 

cosh i9 + cosh (31 



exp 



d^' log (1 + tanh2('"i+i) f ) 
27 cosh(?9 - 



(3.34) 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



The leading term is the energy of mi massive particles with zero momentum and the second 
term describes the finite-size vacuum polarization in the presence of m2 massless particles 
with momenta given by ( |3.37|) . Notice that the contribution from three particle interactions 
in the Yang-Lee model is absent here since the relativistic kink particles do not have 
bound states. The quantum number mi giving the number of massive particles can be 
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identified as the number of domain walls or kinks in the configurations of the RSOS off-critical 
A4 lattice model in Regime III. This identification is possible because the same classification 
of eigenvalues in terms of (m, n) systems applies to the lattice model throughout the off- 
critical Regime III. Indeed, an (m, n) system appears in (2.5b) and (2.6) of |T^ with t = mi, 
s = m2, r = rii and p = n2. Although the periodic case was considered in this previous paper 
the same (m, n) system applies for the (r, s) = (1, 1) boundary condition with double row 
transfer matrices. The identification of the number of massive particles mi with the number 
of domain walls can therefore be easily made by looking at low temperature expansions. 



4 Massive Numerics 

Away from the UV and IR limits, the excited TBA equations and associated auxiliary 
equations can only be solved numerically. While simple iteration of pseudo-energies is usually 
enough for ground-state TBA analysis, extra complications arise for excited TBA equations 
due to the presence of the zeros. 

Our numerical algorithm is to iteratively update the pseudo-energies with previously 
determined values of the zeros and then to find new values for the zeros using the updated 
pseudo-energies. This iteration continues until we obtain the required data with a desired 
accuracy. One delicate point arises when one solves the auxiliary equations. There is no 
natural way to rearrange the second equation in ( p. 271) to express the zeros in strip 2 directly 



„(2) 

in terms of other quantities. Instead, we use the log term ^^'^^ = —i log [tanh^ as 

(2) 

iteration variable which gives, in turn by inversion, the improved values of the zeros in 
strip 2. This ^^'^^ is naturally interpreted as a phase factor. We coded the algorithm in the 
MATLAB™/Octave programming language. Typical running time on a 500 MHz computer 
to achieve an accuracy of five decimal digits is about one minute for a given value of R. 

For the purposes of plotting numerical data it is convenient to normalize the scaling 
energies RE{R). As we have discussed, the leading term in RE{R) diverges linearly with R 
as -R — > 00 while it approaches a constant as i? — 0. To plot the whole flow from UV to IR 
in one plot, we use a normalized scaling function 

^(R) = (4-1) 
^ ^ 2(7r + mi?) ^ ^ 

with the UV and IR limits 

7 

^(0) = lim S{R) = -— + Ar,s + n, S{oo) = lim £{R) = mi. (4.2) 

-R— >0 240 ' R^co 



4.1 Sector (r, s) = (1, 1) 

In Figure 5 we show our numerical results for the (r, s) = (1, 1) boundary condition. This 
sector has an even number of zeros in each of the two strips. The vertical axis is the 
normalized scaling function S{R) and the horizontal axis is log^Q (mi?). We plot selected 
normalized scaling energies for up to mi = 6 zeros in strip 1 and m2 = 2 zeros in strip 2 for 



19 



mi 


m2 






A + n 


£{oo) 


mi 


m2 




# 


A + n 


£:(oo) 











L 1 J 








2 





9 


[3151 

L 1 J 


11 


2 


2 








L 1 J 


2 


2 


4 





3 


[1|31 

L 1 J 


11 


4 


2 





1 


L 1 J 


3 


2 


4 


2 


5 


[1|6] 

L 1 J 


11 


4 


2 





2 


[2121 

L 1 J 


4 


2 


2 





10 


[1161 

L 1 J 


12 


2 


2 





3 


[2121 

L 1 J 


5 


2 


4 





4 


fl|5l 

L 1 J 


12 


4 


2 





4 


[313] 

L 1 J 


6 


2 


4 


2 


6 


[1191 

L 1 J 


12 


4 


4 


2 





L 1 J 


6 


4 


2 





11 


[1161 

L 1 J 


13 


2 


2 





5 


[3131 

L 1 J 


7 


2 


4 





5 


[1161 

L 1 J 


13 


4 


4 


2 


1 


L 1 J 


7 


4 


4 


2 


7 


[1|11] 

L 1 J 


13 


4 


2 





6 


[4|4] 


8 


2 


2 





12 


[1|7] 


14 


2 


4 










8 


4 


4 





6 


[1|9] 


14 


4 


4 


2 


2 


[2|2] 


8 


4 


4 


2 


8 


[1|15] 


14 


4 


2 





7 


[4|4] 


9 


2 


6 


2 







14 


6 


4 





1 




9 


4 


2 





13 


[1|7] 


15 


2 


4 


2 


3 


[3|3] 


9 


4 


4 





7 




15 


4 


2 





8 


[3|5] 


10 


2 


4 


2 


9 


[1|17] 


15 


4 


4 





2 


[2|2] 


10 


4 


6 


2 


1 


[1|2] 


15 


6 


4 


2 


4 


[1|5] 


10 


4 















Table 1: Energy levels for the (r, s) = (1, 1) boundary condition with A = 0. The quan- 
tum numbers mi, m2, E-'^ — E-'^^^'' + E-'^^^^ shown along with the conformal-massive 
dictionary connecting the UV conformal data A + n with the number of massive particles 
S{oo) = rrii in the IR limit. The degeneracies ^ = [^i|^2] indicate that £i levels are plotted 
out of the £2 levels with the given quantum numbers. The plotted energy levels are complete 
in the conformal limit up to n = 9. 
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■1 1 

log mR 



Figure 5: Normalized scaling energies £{R) 
(r, s) = (1, 1) boundary condition. 



2(^+mR) plotted against logj^Q mR for the 
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.y I I I I I I I I 

-1 -0.5 0.5 1 1.5 2 2.5 

log mR 

Figure 6: Locations of zeros plotted against log^Q mi? for a typical energy level for the 
(r,s) = (1, 1) boundary condition. Here rm = 6, ma = 2 and = (1, 0, 0, 0, 0, 0|1, 0). 

As -R — s> oo, the six zeros in strip 1 approach zero while the two zeros in strip 2 converge 
to constants given by ( p.37] ). Since the zeros are in different strips there are no collisions of 



allowed quantum numbers including the lowest 30 levels. Note that eigenvalues which are 
degenerate in the conformal UV limit become non-degenerate when the perturbing thermal 
field is turned on and the conformal symmetry is broken. Table 1 summarizes how the UV 
descendant levels flow into the IR particle states. The plotted energy levels are complete 
in the conformal limit up to n = 9 corresponding to the expansion of the finitized Virasoro 
character ( |2.4(J| ) in the limit N ^ oo 



oo 
2 



oo 
4 



oo 
4 



+ q 
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OO 

6 



1 + + g + + 2g=^ + 3q^ + + 4g6 + 4g^ + . . . ) 

+ g^(l + q + ...) + q\l + g + 2^^ + 3g^ + ...) + . . 
l + q^ + q^ + 2q^ + 2q^ + 4g^ + 4g^ + 7g^ + 8g^ + . . . 



+ q 



14 



oo 
6 



+ . . . 



(4.3) 



The location of zeros also flow as R changes. To illustrate this, we plot in Figure 6 
the locations of the eight zeros for the energy level with mi = 6, m2 = 2 and quantum 
numbers = (1, 0, 0, 0, 0, 0|1, 0). Near the UV limit, the zeros are linear in log{mR) 

as expected in (|3.3CI| ). One can check that the constant separations between them agree 
with those from the "conformal TBA" equations. The IR behaviour of the zeros is also 
as expected. The six zeros in strip 1 decay exponentially to as plotted against log(mi?). 
This is in accord with the fact that the zeros approach as 1/mR and the constants of 
proportionality are conflrmed to be precisely those in ( |3.35| ). The massive particle states 
become frozen in the IR limit and accordingly the two zeros in strip 2 converge exponentially 
to the flnite constant values given by ( |3.37| ). The massless particles have prescribed momenta 
in the same limit. These behaviours are reproduced consistently in the analysis of all other 
boundary conditions and quantum numbers. 

4.2 Sector (r, s) = (3, 1) 

In this sector, the (m, n) system is 



m + n = - [{N + l)ei + Im] 



(4.4) 



and N and mi are even while m2 is odd. Hence there is an even number of zeros in strip 1 
and an odd number in strip 2. Repeating the derivation of the TBA equations leads to the 
same equations as for the (r, s) = (1,1) boundary condition but with e"''^^''^ replaced with 
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-e fjj^g TBA equations in this sector thus become 



mi 



(1) 



1 

2^ 



log tanh - - 2^ log ^tanh ( 2 + 2 ^ ^ ^2 
, log(l - e-^^^'^')) 



(1) 



dif 

cosh(i? - 1?') 



§ ^ r § ^ ^(2) 

e2(^?) = 2mi?cosh ^? - log tanh^ zJ tanh (- + tanh (- - -|-) 

fc=i 



,log(l + e2^ 
cosh(T? - ■&') 

The auxiliary equations are 



1 /"^ 
27r 



2mi? sinh 



m2 

— i log tanh ^ 



« '"8 _ , g log [ ta.h ( ^ + 

fc=l 



27r sinh(/3f ^ - ^9) 



)] 



it=i 



/3f ) + /?) 



-)J -ilog 



0= :- 



4 2 

rf?? log(l - e-'^{d)) 



2 /TT^ 



tanh^ ^ 



mi 



27r smk{l3f^-d) 



-i J]loj 



tanh (- + ^ 



- i ^ log 



tanh 



(1). 



4 2 
and the scaling energy becomes 



— i log 



tanh^ f ^ - ^ 
V 4 2 



(2) 



)] 



I (2) 



mi 



= 2mR ^ cosh /?] 



(1) 



mR 



/oo 
(it? cosh "Q log (1 — e' 
■00 



(4.5) 



(4.6) 



(4.7) 



In Figure 7 we show our numerical results for the (r, s) = (3, 1) boundary condition. The 
vertical axis is the normalized scaling function S{R) and the horizontal axis is logj^gl'^-^)- 
We plot selected normalized scaling energies for up to mi = 4 zeros in strip 1 and m2 = 1 zero 
in strip 2 for allowed quantum numbers including the lowest 49 levels. Table 2 summarizes 
how the UV descendant levels flow into the IR particle states. The plotted energy levels arc 
complete in the conformal limit up to n = 9 corresponding to the expansion of the finitized 
Virasoro character in the limit N ^ 00 
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Table 2: Energy levels for the (r, s) = (3,1) boundary condition with A = 3/2. The 
quantum numbers mi, m,2, E-^ = E-^*^^'* + E-^''^'' shown along with the conformal- 
massive dictionary connecting the UV conformal data A + n with the number of massive 
particles S{oo) = mi in the IR limit. The degeneracies 7^ = [ii\i2] indicate that ii levels are 
plotted out of the £2 levels with the given quantum numbers. The plotted energy levels are 
complete in the conformal limit up to n = 9. 
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Table 3: Energy levels for the (r. ,s) = (2,1) boundary condition with A = 7/16. The 
quantum numbers mi, m2, E ~ E + E-^'^^'' shown along with the conformal- 
massive dictionary connecting the UV conformal data A + n with the number of massive 
particles £(00) = mi in the IR limit. The degeneracies # = [^i|^2] indicate that ii levels are 
plotted out of the £2 levels with the given quantum numbers. The plotted energy levels are 
complete in the conformal limit up to n = 6. 
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Figure 7: Normalized scaling energies S{R) = 2(^+mR) plotted against log j^q mi? for the 
(r, s) = (3, 1) boundary condition. 
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4.3 Sector (r, s) = (2, 1) 

In this sector, the (m, n) system is 



m + n = - {Nci + 62 + Xm) 
2 



(4.9) 



and N, mi and 777,2 are all odd. Hence there is an odd number of zeros in strip 1 and in strip 2. 

Repeating the derivation of the TBA equations leads to essentially the same equations as for 



the (r, s) = (1, 1) boundary condition but with e ^^^"^^ replaced with — e" 



1 , 2 and an 



additional tt in the second auxiliary equation. The TBA equations in this sector are in fact 

^ r . 



ei(^) 



— log tanh^ y. log tanh ( — I — tanh ( 

27r cos] 



3h(i? - 1?') 



62(1?) = 2mi?coshi? - logtanh^ - - log | tanh (- + -|-) tanh (- - -|-) 



(4.10) 
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cosh(i? - 

and the auxiliary equations 



2mi?sinh/?j^^ 



y 2^ 



di? l og(l-e-^i(7^)) 
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di? log(l-e-^2(^)) 
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4 2 



(4.11) 



-z^log 



V 4 2 



)]-^log[tanh^ (t - 2 )] + ^^^'^ + 



with the scaling energy 



'"I 77?/? /""^ 

i?E(i?) = 2mi? V cosh /5f ' / cosh log (1 - 



-£2(1?) ^ 



(4.12) 



In Figure 8 we show our numerical results for the (r, s) = (2, 1) boundary condition. The 
vertical axis is the normalized scaling function S{R) and the horizontal axis is log^ol^^^-^)- We 
plot selected normalized scaling energies for up to mi = 5 zeros in strip 1 and m2 = 3 zeros 
in strip 2 for allowed quantum numbers including the lowest 18 levels. Table 3 summarizes 
how the UV descendant levels flow into the IR particle states. The plotted energy levels are 
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Figure 8: Normalized scaling energies S{R) = 2(^+mR) plotted against log j^q mi? for the 
(r, s) = (2, 1) boundary condition. 
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complete in the conformal limit up to n 
Virasoro character in the limit N ^ oo 



6 corresponding to the expansion of the finitized 
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+ ... 

{1 + q + + + + + + q'^ + + + q^^ + . . .) 
+ + q + 2q^ + 3q^ + V + 5q^ + 7q^ + 8q'^ + . . . + q) 
+ q\l + q + ...) + g^°(l + ...) + ••■ (4.13) 

1 + g + + 2g^ + 3g^ + V + + 8/ + lOq^ + Uq^ + 18g^° + . . . 



5 Discussion 



In this paper we have derived and solved numerically the TBA equations for all excitations 
for the massive tricritical Ising model with boundary conditions labelled by (r, s) = (1, 1), 
(2, 1) and (3, 1). The analysis can be extended to the other primary boundary conditions 
(r, s) = (2,1), (2,2) and (3,2) by allowing for frozen zeros and introducing two (m, n) 
systems in the classifications of eigenvalues. It would also be of interest to extend our 
analysis to periodic boundary conditions. The main new feature of such a calculation would 
be the classification of the periodic eigenvalues which would entail many (m, n) systems and 
must allow for different patterns of zeros in the upper and lower half planes related to the 
two (left and right) copies of the Virasoro algebra. If our analysis was extended to periodic 
boundary conditions it would allow a direct comparison of the results of our lattice approach 
with the results of the Truncated Conformal Space Approximation (TCSA) which are good 
for small niR. 

In Part II of this series of papers we will derive and solve numerically the TBA equations 
for all excitations for the massless flow from the tricritical to critical Ising model. This has 
some interesting additional features because zeros can collide during the flow leading to 
changes in the classification of eigenvalues and to a flow between Virasoro characters of the 
two theories. 
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